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Abstract 

The problem of finding most general form of the classical integrable relativistic 
models of many-body interaction of the BCn type is considered. In the simplest 
nontrivial case of n=2, the extra integral of motion is presented in explicit form 
within the ansatz similar to the nonrelativistic Calogero-Moser models. The result- 
ing Hamiltonian has been found by solving the set of two functional equations. 
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The one-parameter extensions of the Calogero-Moser integrable particle systems first 
proposed by Ruijscnaars and Schneider [1,2] have been extensively studied during last 
decade [3,4,6-17]. It turned out that both classical and quantum versions of these sys- 
tems (often called "relativistic" but, as explained in [6], this term is misleading) were 
connected with various branches of mathematical physics. The quantum RS models with 
trigonometric type of interaction have been found to be relevant in the theory of Mac- 
donald polynomials and the integrability of their elliptic version has been investigated 
within the Yang-Baxter approach [8,9]. The apphcations of classical models range from 
characteristics of multisoliton solutions of the sine-Gordon equation to the Seiberg-Witten 
theory [13]. 

The original n-particle RS model has the Hamiltonian 

n n 
j=l k^j 

where p,q are canonical variables with Poisson brackets {pj,Pfe}p = {?j,?fe}p = 0, 
{Q.j-iPk\p — ^jk- It is translation invariant and can be associated with An-i root sys- 
tem. Extensions to other root systems were proposed [7,10-17] in ways similar to the 
Olshanetsky-Perelomov construction for Calogero-Moser models [18]. One can find the 

first of them in the paper by Schneider [4] who considered the non-translation-invariant 
model with an external field and the Hamiltonian 

n 

Hs^Hn + Y.W{q,). 

The Schneider model has been found to be integrable for hyperbolic choice of the two- 
particle interaction, 

F2(g) = 1 + /sh2(ag), W{q) = eycp{±2aq) , 

which coincides upon the limit 

Pj Pj/c, Qj QjC, g g/c, A A/c, c^oo 

with well-known Adler system [5]. The next step has been done in the papers by van 
Diejen [7] who has studied the quantum RS systems and has found their nine-parametric 
extension (with one constraint) of the BCn type in the most general elliptic case. The 
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results of [7] have been proved rigorously for n = 2 only. Komori and Hikami [10,11] 
also considered the elliptic potentials starting from the paper by Hascgawa [9] within 
the framework of quantum inverse scattering method and gave a rigorous proof of the 
commutativity of van Diejen operators for arbitrary n and with no constraint among 
couplings. As for classical models, the Lax operators were obtained in explicit form only 
for very particular case of Z2 folding of the original and ^2n+i RS model [14-16]. 

The aim of this Letter is to find most general form of classical integrable extensions 
of the RS hamiltonian of the SC-type. The general procedure of constructing the Lax 
pairs for these systems is still lacking, so we shall consider the simplest but nontrivial 
systems of two degrees of freedom. In this two-particle case, one needs to find only one 
extra integral of motion J{pi,P2, Qi, 12) so as to prove the integrability of the system. 

We shall consider the following Ansatz for the two-particle Hamiltonian and an extra 
integral of motion J, 

H2 = Vi{qi, q2) chpi + V2{qu ^2) chps + 92), (1) 

J = A(gi, gs) chpi chp2+n{qu ^2) shpi shp2+Pi(gi, ^2) chpi-Fp2(gi, ^2) chp2+'r(gi, ^2), (2) 

where Vi_2, W , A, pi^2 and r are yet unspecified functions of particle coordinates. We 
shall suppose that H2 and J are symmetric with respect to the permutation of particles, 

V2(gi, 52) = Vi{q2, qi), W{qi, ^2) = W{q2, qi), (3) 

52) = A(g2, gi), ^2) = /x(g2, gi), P2(gi, ^2) = Pi(g2, gi), r{qi, ^2) = r(g2, qi)- 

(4) 

The calculation of the Poisson bracket {H2,J}p and equating it to zero results in an 

overdetermined system of 10 equations for 8 unknown functions, 

dX dV, dV^ dX dV2 .dV2 ^ 

Vi- A- fi— = 0, V2- /i- A— = 0, (5) 

oqi oqi oq2 oq2 oqi oq2 

dii dVi dVi dii dV2 dV2 _ 

Vi- p.- A— = 0, V2^ A- /i— = 0, (6 

oqi oqi oq2 oq2 oqi oq2 

dpi dVi dp2 dV2 , . 

Vi- pi— = 0, V2^ p2^ = 0, 7) 

oqi dqi dq2 dq2 

dr da dW dr da dW , , 

Vi—-V2^-Pl^ = 0, V^2T^-V^17T^-P2T^ = 0, (8) 



dqi dq2 dqi dq2 dqi dq2 

dp2 dV2 dW dW ^ ^dpi dVi dW ^dW 



To solve the system (5-9), it is useful to start from four equations (5,6). They can be 
easily linearized and the general solution which takes into account also the symmetry (3) 
reads 



= F{q^ - q2)G{q^ + q2)R{qi), = F{q, - q2)G{q, + q2)R{q2), (10) 



where F, G, R are (still unknown) functions of one argument, and 

G\q, + q2)]R{qi)R{i 

11) 



A = -hF'(q,-q,)+G\q, + q2)]Riq,)Riq2), /x = l[F\q,-q2)-G'iq, + q2)]R{qi)Riq2). 
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The solution of (0) has the form 

Pi = Vid{q2), P2 = V2d{qi). (12) 
The equations (H) can be rewritten as 

These equations with the use of (|TTD can be cast into the form 

^ = R{q2)[lR'{q2){F' - G^) - R{q2){FF' + GG')] + ^(^2)^, (M) 

^ = R{q,)[\R'{q,){F' - G') + R{q{){FF' - GG')] + ^(^1)^, (15) 

where prime means differentiation of the function with respect to its argument. Now the 
equality of second mixed derivatives of r obtained from (0,|l^) yields 

— — -(rf(g2) - d{q{)) + d\q2)^ - d'{qi)—- 
oqiOq2 oqi dq2 

+ ^[\R{q2)R\q2){F' - G') - R\q2){FF' + GG')] 

^ []-R{q,)R!{qi){F^ - G^) + R\qi){FF' - GG')] = 0. (16) 



dqi'2 

It is useful to introduce instead of F, G and R new functions by 



:l+/)l/^ G = (l+^?)l/^ R={i+r)'/'. (i?) 



The substitution of ( p!0| - |T^ , p!7D into the last pair of equations (|^) yields after some trans- 
formations 
dW 1 

^ = -[{f + 9'){d{q2)-d{q,)) - (1 + f){d'{q,)+d'{q2)) + {l+9){d'{q2) - d'{q,))], (18) 



dW 1 

^ = 2 [(/' - 9'){d{q2) - d{q,)) - (1 + f){d'{q,) + d\q,)) - (1 + g){d\q,) - d\q,))]. (19) 
The compatibility condition for the equations (p!8| - p!9D reads 

2(/"(gi - q2) - 9"{qi + q2)){d{q2) - d{q,)) + {f{q, - q^) - g{q, + q2)W {q^) - d"{q,)) 
-3/'(gi - q2){d'{qi) + d\q,)) + ^{q, + q2){d\q^) - d\q^)) = 0. (20) 

The equation (p!6| ) can be transformed with the use of (p!7| - p!9D to the form 

2(/"(gi - q2) - 9"{qi + q2)){r{q2) - r{qi)) + (/(gi - ^2) - g{qi + q2)){r"{q2) - r"{qi)) 
-3/'(gi - g2)(r'(gi) + r (ga)) + 3g'{qi + q2){r'{qi) - /(ga)) 
= (rf(gi) - d{q2))[{d"{q,) + d"{q2)){g{qi + ga) - /(gi - g2)) + 3/'(gi - g2)(rf'(g2) - rf'(gi)) 
+3^7'(gi + g2)(c^'(gi) + c^'(g2))] + 2{g{q^ + ga) - /(gi - g2))(rf''(gi) - ^^''(g2)). (21) 

One is left with two functional equations (^, ^1]) . The first of them allows one to find the 
functions /, g and d while the second one can be regarded as inhomogeneous functional 
equation for the function r(g). The most general solution to (^) with f = g has been 
found in [19], 

f{q)=9{q) = fop{q), (22) 
where p{q) is the Weierstrass function with two arbitrary periods 2lji,2lj2, 

3 

d{q) = Y.(^ip(l + ^i)^ (23) 

i=0 

where we use the notation loq = 0, = —ui — uj2, and {aj} are four arbitrary real 
parameters. 

This solution has to be used in equation ( pi]) which is bilinear in d{q). Plugging ( |2^ 
P5D into right-hand side of (^Tj) results, besides the terms quadratic in {a,}, in 6 terms 
proportional to aia^ with i ^ k. However, it is possible with the use of Mathematica 
to show that all these cross-terms vanish and the right-hand side of ( PT| ) acquires simple 
form 

3 

8lI«»V(9i +^i)p'(92 + ^i)b(g2 + c^i) - p{qi+uJi)]. 

1=0 

It allows one to find the general solution to (|2l|), 

3 

r{q) = T.Mq + ^i) + (y-p\q + ^i)], (24) 

i=Q 



where {rj} is yet another set of four arbitrary real parameters. After solving equations 
(|T^-|T^) the resulting Hamiltonian reads 

H, = j2 chp^ n{[l + fi^J - + fiQ, + QkW + r{qj)]V^' + W{q,, g^), (25) 

j=i k^j 



where /, g and r are given by and 



W{qi, q2) = - Y.[fo<^ip(qi + ^i)pi.q2 + i^i)] - d{qi) - d{q2). 



(26) 



j=0 



In the papers [7,10,11] the Hamiltonian was written in terms of the products of the 
Weierstrass sigma functions. Our solution (P^ ) for 1 + r(g) can be also cast into this form 
by using the formulas 



p{q)p{q + ^i) 
p{q) - p(a) 



p{uJi)[p{q) + p{q + Ui)] - p'^iuJi) + p"{uJi)/2, 
cr(a — q)a{a + q) 



a'^{q)a'^{a) 

The second invariant J is determined by the formulas (|lll-0) and (0-|l5|). The explicit 



expression for the function r(gi, ^2) can be also obtained by the integration of (14-15). In 
the elliptic case, it reads 



r{qi,q2] 



-^^fiqi - q2) + /(gi + ^2)] - ^(^1)^(^2) 

3 



i=0 



n pi.qi + ^i)pi.qi + ^k), 



where /(g) and d{q) are given by (^2]-^) and 

piqi)p'^iq2) + piq2)p'\qi) 



(27) 



^(91,92) = p{qi)piq2 

X(?l,?2) = 



4(p(gi) - Piq2)y 

pHqi)p'Hq2) + pHq2)p'Hqi) 



4(p(gi) - p{q2W 

Taking correct hyperbolic (trigonometric) and rational limits of (22-26) is a nontrivial 
task. To do this, we have solved again the equations ([2T| - |22D with the use of Mathematica. 
The results are given by (^Sj) with the following sets of functions f,g,d,r,W and r: 



(a) Hyperbolic case: 

/(g) 
d{q) 

r{q) 



/o 

+ 



^0 ^ n 
sh^aq ch^aq 



+ a2 ch2ag + ^3 ch4ag, 
+ r2 ch2ag + ch4ag 



+ 



22 2 

H — — h Q;2a3 cnoag H — ^ cnsag, 



sh aq ch ag 



"0 



sh^agi sh^ag2 ch^agi ch^ag2 



+ 4«3 ch2agi ch2ag2 



(28) 
(29) 

(30) 

)■ 

(31) 



-d{qi)d{q2) - ( sha(gi - ^2) sha(gi + ga)) [( ch2agi ch2ag2 - l)/2 
+ro( ch^agi + ch^ag2) + ri( sh^agi + sh^ag2) + r2( sh^agi ch2agi ch^ag2 
+ ch^agi sh^ag2 ch2ag2) + r3( sh^2agi sh^2ag2 + sh^agi ch^ag2 + sh^ag2 ch^agi) 
+aQ(coth^ agi + coth^ ag2) + a^(tanh^ aqi + tanh^ ag2) + a2a3( ch2agi + ch2ag2) 
x(4ch2agi ch2ag2 - 3)(ch2agi ch2ag2 - l)/4 + a3(8ch^2agi ch^2ag2 
+9 ch2agi ch2ag2 — 4 ch2agi ch2ag2( ch2agi + ch2ag2)^ — l)/4 
— /o 16a;o«i( sh2agi sh2ag2)~^ + aoa2( shagi shag2)^^(l + 2sh^agi + 2sh^ag2) 



+aia2{ chagi chag2) ^(1 — 2 ch^agi — 2 ch^ag2) + a;o«3(16( sh^agi + sh^ag2) 



+8(ch agi sh ag2 + ch ag2 sh aqi) + (shagishag2) ) +aia3(16(ch aqi 
+ ch^ag2) — 8( sh^agi ch~^ag2 + sh^ag2 ch~^agi) — ( chagi chag2)~^) 



(32) 



where {r,, ctj} is the set of 8 arbitrary real parameters which does not coincide with those 
introduced for general elliptic case in (p3D , (l2^) (the same remark holds also for subsequent 
rational and asymmetric hyperbolic models), 
(b) Rational case: 



d{q) 
r{q) 



9{q) = foq~ 



aoq ^ + aiq^ + Q;2g'^ + a3g'', 
f'oq''^ + riq^ + r2g^ + r3g^ 



, 2„12 



+aog + (2a;ia3 + Q;2)g + 2a2Q;3g + "3? 



(33) 
(34) 

(35) 
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W{qi,q2) = -^d{qi) - fo[ao{qiq2) + a2{ql + ql) + a3{qf + q2 + 3qlql)], (36) 



i=l 

-riqi, q2) = [f{qi - ^2) + /(gi + ^2)] - d{qi)d{q2) 



2 







o lO 22, 22,0 44. /2| 2\/ 2 -2 -2 I 2 4 4 



+Oilqlql + a2a3qfq^{ql + q^)) + aiasqfqUqf + ^2 + gig2(gi + ^2)) 
-/o [aoai(gr^ + g^^) + aoa2(gig2'^ + gki"^) 

+aoa3(gig2"^ + g2gr^ + 3(gi + iD) + aia2qfql\ , (37) 

where /(g) and d{q) are defined by (p3|-p^) and {aj,rj} is the set of 8 real parameters 
which do not coincide with those in (0). There is yet another hyperbohc degeneration 
which corresponds to the generalization of Schneider case [4] of interaction with external 
field, 

/(g) = -r^, 9 = 0, (38) 
sh aq 

d{q) = aoe-^^^^ + aie^"" + a2e-'«9 + aae'"^ (39) 
r(g) = + rie^'^" + rse-^'^" + rge^'^'? 

+2aoa2e-^"^ + 2aia3e^"^ + a^e"®"'? + ale^''\ (40) 

W^(gi,g2) = -E%i)-2/o(«2e-'"^«^+''^) + «3e'"('^^+«^)), (41) 

i=l 

r(gi,g2) = -%i)rf(g2) -/o[(2sh2a(gi -g2))-'[l + ro(e-2'^^^ +e-2'^'?^)/2 

+ri(e^''^^ + e^'"^^)/2 + r2e-2"(9i+92) + ^ggSata+ga) ^ aga2e~^"^'^^+''^)(e"2'''?^ + e"^"''^) 

+2«o«3(e''^''i + e^"""') + 2ai«2(e-'"^^ + e'^^^^) + 4«2«3 ch2a(gi - g2)].(42) 
The original Schneider's result [4] for two particles corresponds to the case of all rj=0, 



^2 = 03 = in (IH-ll 



Let us summarize our results. We have demonstrated the Liouville integrability of the 
classical Ruijsenaars- Schneider Hamiltonian of the BC2 type. It contains nine arbitrary 
coupling constants with no constraints in the general elliptic case and its hyperbolic and 
rational degenerations. Its form (^) is suitable for taking "nonrelativistic" limit p p/c, 
q cq, c ^ 00: it is necessary to rescale /o fo/c^, fi — ^ rj/c^, ctj for getting most 
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general five-parametric two-particle Calogero-Moser Hamiltonian given in [19]. It would 
be of interest to construct the Lax representation for our general nine-parametric case so 
as to prove the classical integrability for arbitrary number of particles. At present time 
it was done only for very particular cases [14-17]. 
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tially supported by the Grant-in aid from the Ministry of Education, Culture, Sports, 
Science and Technology, Japan, priority area (#707) "Supersymmetry and unified theory 
of elementary particles" . 

References 

1. S.N.M. Ruijsenaars and H. Schneider. Ann. Phys. (NY) 170 (1986) 370. 

2. S.N.M. Ruijsenaars. Commun. Math. Phys. 110 (1987) 191; 115 (1988) 127. 

3. M. Bruschi and F. Calogero. Commun. Math. Phys. 109 (1987) 481. 

4. H. Schneider. Physica D26, (1987) 203. 

5. M. Adler. Commun. Math. Phys. 55 (1977) 195. 

6. H.W. Braden and R. Sasaki. Prog. Theor. Phys. 97 (1997) 1003. 

7. J.F. van Diejen. J. Math. Phys. 35 (1994) 2983; 36 (1995) 1299. 

8. A. Antonov, K. Hasegawa and A. Zabrodin. Nucl. Phys. B503 (1997) 747. 

9. K. Hasegawa. Commun. Math. Phys. 187 (1997) 289. 

10. Y. Komori and K. Hikami. J. Phys. A30 (1997) 434. 

11. Y. Komori and K. Hikami. J. Math. Phys. 39 (1998) 6175. 

12. K. Hasegawa, T. Ikeda and T. Kikuchi. J. Math. Phys. 40 (1999) 4549. 

13. H.W. Braden, A. Marshakov, A. Mironov, A. Morozov. Nucl. Phys. B558 (1999) 
371. 

14. B. Hou, W.-L. Yang. J. Math. Phys. 41 (2000) 357. 



9 



15. K. Chen, B. Hou and W.-L. Yang. J. Math. Phys. 41 (2000) 8132. 

16. K. Chen and B. Hou. The D„ Ruijsenaars-Schneider Model, |hep-th/010203(: 



17. J. Avan and G. Rollet. Structures in BCn Ruijsenaars-Schneider Models, 



nlin.SI/0106015 



18. M.A. Olshanetsky and A.M. Perelomov. Phys. Reports 71 (1981) 313. 

19. V.I. Inozemtsev. J. Phys. A17 (1984) 815. 



10 



